On the valence bond solid in the presence of Dzyaloshinskii-Moriya interaction 
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We examine the stability of the valence bond solid (VBS) phase against the Dzyaloshinskii-Moriya 
(DM) interaction in the bipartite lattice. Despite the VBS is vulnerable against the antiferromag- 
netic interaction, for example in the Q-J model proposed by Sandvik, where the quantum phase 
transition occurs at J* /Q — 0.04, we found that on the contrary the VBS is very stable against the 
DM interaction. The quantum phase transition does not occur until D/Q goes to infinity, where D 
is the strength of the DM interaction. The VBS in the ALKT model and the Haldane gap system 
also exhibit the same property. 
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Valence bond solid (VBS) is the array of singlets that 
all the antifciTomagnctic spins are paired statically with 
one of their nearest-neighbor spins in the real space. It 
preserves the spin rotational symmetry because all spins 
are paired in singlets and consequently the total mag- 
netization is zero. Since the spins do not change their 
pairing partners , those singlets usually form a periodic 
array that breaks the translational symmetry in addi- 
tion to the lattice symmetry. Its "liquid phase", usu ally 
recognized as the resonating valence bond (RVB)^, re- 
stores the translational symmetry by forming the dynam- 
ical singlets that spins can change their pairing partners 
dynamically with the infinite or the short range pairing 
strength. Both of them have attracted great attention to 
the condensed matter physicists because of their unusual 
ground state and the excitation properties. VBS usu- 
ally accompanies with a spin gap by breaking one pair 
of singlets. The RVB, on the other hand, often requires 
the critical spin correlation, and interestingly the fun- 
damental excitations may carry the fractional quantum 
number. 

Albeit they are immensely interesting, pure spin mod- 
els to host them as the ground state is not easy to con- 
struct, especially in the dimensions greater than one. 
Therefore, their ground state as well as the excitation 
properties are not yet well understood. Part of the rea- 
sons is that the antiferromagnetic energy often dominates 
when competing with the VBS and the RVB states. Con- 
sequently, the ground state favors the antiferromagnetic 
(AFM) state rather than the VBS or the RVB states. 

Since we are interested in the properties of the VBS, 
the model Hamiltonians to host it as the ground state 
are necessary. The easiest way to realize the VBS ground 
state is the one-dimensional spin Peierls system, where a 
soft lattice is required as well as the strong spin-lattice 
interaction. The spin Peierls system does not serve our 
satisfaction, because we are looking for the VBS state 
from a pure spin model with low spin quantum number 
at best. Another possible system with the VBS ground 
state in one dimension is the Haldane gap system, namely 
the Heisenberg model with integer spins. One realization 
of the Haldane gap system is the AKLT modeP described 



by the following Hamiltonian 

Haklt = ■ $+1 + ■ S t+1 ) 2 + |], (1) 

i 

with S = 1 and J > 0. The spin one on each site is 
constructed by the symmetric projection of two spin-1/2 
spins. Since the coordination number is two, these two 
spin-1/2 spins on every site form the singlets with each 
nearest neighbor. This state is certainly the VBS state 
and can be proved to be the exact ground state of Eq.Q. 

In two dimensions, there are few spin models that 
exhibit VBS ground states. For example, the SU(N) 
Heis enber g model is shown to have the columnar VBS 
stattP^ for N > 5 and the AFM for N < 5. Recently, 
Sandvik^ proposed a ring-exchange-like model and con- 
sidered the quantum phase transition between the VBS 
and the AFM, described by the following Hamiltonian 

h qj = -q (s l -s J - 1 -)0 k -s l - 1 -) + jJ2s l -s J (2) 

<ijkl> <ij> 

where < ij > denote the nearest neighbor spins and 
< ijkl > are the coordinates of a plaquette with the 
counterclockwise or the clockwise order. We shall call it 
QJ model in the rest of the paper. Using this model, 
Sandvik found that a second order quantum phase tran- 
sition from the VBS state to the AFM state take place at 
J/Q = 0.04. This small ratio of J/Q indicates that the 
VBS state is indeed vulnerable to the AFM perturbation. 
Another model that might host the VBS ground state is 
the Heisenberg model in the kagome lattice^. Numeri- 
cal calculation showed that it is a VBS with 36-site unit 
cell. Although the final verdict of the ground state of the 
Heisenberg model in the kagome lattice is still under de- 
bate, this model is worth mentioning, because it relates 
to the ZnCu 3 (OH) 6 Cl2 that shows no spin ordering down 
to tens of the mK as well as the absence of the spin gap 
down to 0.1me\^. The lack of the spin ordering is also 
one of the important properties of the VBS state. How- 
ever, since there is a spin gap in the VBS system, the 
Heisenberg model can not fully describe the material. 

In this paper, we study the Dzyaloshinskii-Moriya 
(DM) perturbation in the VBS state. Part of the motiva- 
tion was enlightened by the experimental fact that VBS 
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systems often occur in the lattice with low crystal sym- 
metry for example the Ni^HgNa^NC^ClCU) (NENP) 
for the Haldane gap system and the ZnCu 3 (OH) 6 Cl 2 . Es- 
pecially, if the inversion symmetry is absence, the DM 
interaction is often inevitable. Another motivation is 
tempted by the study of ZnCu3(OH) 6 Cl2. Further nu- 
merical study had considered the addition of the DM 
perturbation to the Heisenberg model in the kagome lat- 
tice. Because the kagome lattice does not have the in- 
version center, the DM interaction must play some role 
in this material. ESR measurement detect the strength 
of the DM interaction is roughly 0.08J^2J On the other 
hand, theoretical study-^ suggests that there is a quan- 
tum phase transition from the VBS state to the AFM 
state at D/J = 0.1. If it is true, the VBS state would 
be also as vulnerable against the DM perturbation as 
well as the AFM one. Regardless the inconsistency of 
their results with previous studies that strikes our un- 
derstanding of ZnCu3(OH) 6 Cl2, the question "does the 
vulnerability of the VBS by the DM interaction also hap- 
pen in other VBS systems?" has its own weight, which 
will be addressed in this paper. 

To get the general understanding to that question, we 
consider the VBS states in the AKLT and the QD mod- 
els, where the QD model is meant by letting J = and 
adding the DM term in the QJ model. These two mod- 
els are chosen, because they are of the low-spin models 
and there is almost no doubt that the ground state in 
these two models is the VBS state. Our results, on the 
contrary, show the stability of the VBS states against the 
DM interaction in those two models. The quantum phase 
transition does not occur at any finite D. We computed 
the spin gap and found that it does not close until D/Q 
in the QD model (or D/J in the AKLT model) goes to 
infinity, where D is the strength of the DM interaction. 

Our results are consistent with the previous study of 
the XXZ Heisenberg model of spin one. Consider the 
Haldane gap system with the DM interaction 

H=jJ2Si-S i+1 + J2Di-(SiXS i+1 ). (3) 

i i 

If we choose D || z, Eq.|3]) becomes 

H = J2iJ^Si +1 + h.c. + JS?S? +1 ], (4) 



where J± = J + iDi. It implies that the DM term can 
be gauged away and maps to the XXZ Heisenberg model 

Hxxz = £(Sf Sf+i + SfSf +1 + \S*S? +1 ). (5) 



with A < 1. Obviously, the SO (3) symmetry is explic- 
itly broken to U(l) for a finite D. In this case, one 
might expect that the quantum phase transition to the 
XY- phase might occur at some A*. However, numeri- 
cal study^ showed that the quantum phase transition, 
which is shown to be of the Kosterlitz-Thouless type, 



does not happen until A = 0. In other words, only in- 
finitely strong DM interaction can close the spin gap and 
destroy the VBS state. Therefore, very unlike the vulner- 
ability of the questionable VBS state in the Heisenberg 
model in the kagome lattice, the possible VBS state in 
the Haldane gap system is very stable against the DM 
interaction. The key lies in the coincidence that the DM 
term can be gauged away and the spin gap is robust at 
any finite A < 1 in the XXZ Heisenberg model. 

Since the AKLT model was proposed to realized the 
Haldane gap, we expect that the same results will hap- 
pen, and indeed it is. We also examine the VBS state 
in the QD model and find the same stability. One might 
speculate that the vulnerability of the VBS in the Heisen- 
berg model in the kagome lattice might be due to the 
incapability of the DM interaction to be gauged away. 
In the AKLT and the QD model, the DM interaction 
can not be gauged away either. If the VBS state in the 
Heisenberg model in the kagome lattice is proven to be 
true by more detail examination, our results may imply 
that the VBS state in that case is very different from the 
ones we study here. Therefore, this exercise might shed 
light on the characterization of the VBS states. 

A simple-minded explanation of these results can be 
understood as the following. Consider the case that 
there are only two spins with the interaction J Si ■ S% 
and J > 0. The ground state of this case is the sin- 
glet state with the energy —3 J/4, and the excited state 
is the triplet with the energy J/4. Introducing the DM 
interaction of D \\ z rearranges the eigenvalues to be 
(J/4, J/4, l/4(-J-2\AD 2 + J 2 ), l/4(— J+2V-D 2 + J 2 )). 
Two unaffected states are | Tt> an d | II > in the triplet. 
It is obvious to see that introducing any finite D does not 
cause level crossing. In other words, the spin gap does 
not close at any finite D valued 

We organize this paper in the following way. In the 
section [I] we applied the DMRG method to compute the 
spin gap as the function of D/ J in the AKLT plus DM 
interaction. In section [IT] we use the Lanczos method to 
compute the spin gap in the QD model mentioned above. 



Finally, we conclude our results in the section III 



I. AKLT MODEL 

In this section, we compute the spin gap of the AKLT 
model plus the DM interaction given by 



KLT+DAI 



= H A 



LKT + Hdm 



Hdm — ^ Di ■ (Si x Si+i). 



(6) 



(7) 



We consider the uniform DM vector D \\ z here, where z 
is perpendicular to the chain direction. The result of the 
stagger DM vector will be the same. As mentioned above, 
the nearest-neighbored DM interaction can be gauged 
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away in the Heisenberg model, but it is not true in the 
AKLT model because of the (Si ■ Si+i) 2 term. If we con- 
sider only the Eq.([7]) at the classical level, an magnetic 
ordered state with the 4a periodicity will be favored as 
the ground state, where a is the lattice constant. Similar 
to the QJ model, the quantum phase transition from the 
VBS state to the DM ground state might be expected. 
To resolve this, we apply the DMRG method to compute 
the spin gap with the m to be used around 400. We solve 
the problem using the open boundary condition, and the 
system size goes up to 300 spins. 
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FIG. 1: (Color online) The finite-sized scaling of the spin gap 
in the AKLT model. The energy unit is J, which is taken to 
be one in the calculation, (a) D from to 1.6. (b) D from 
1.8 to 2.8. The dash line for D = 2.8 is the guide to eyes. 

In Fig.Q, we show the finite-sized scaling of the spin 
gap for different D/J. Since the excitation of the VBS 
state is to break a singlet, we focus on the triplet exci- 
tations. We extrapolate the scaling and obtain spin gap 
Aqo in the thermodynamic limit for different D/J. The 
result is shown in the Fig.Q. As our DMRG accuracy is 
around 1 x 1.1 x 10 -3 , we do not observe the spin gap to 
close up to D/J ~ 3. It can be seen that the gap drops 
rapidly around D/J ~ 0.6 and is followed by a slow sat- 



uration, in which the gap remains finite and shows no 
evidence of the quantum phase transition at finite D/J. 
Although it is not shown here, we have checked that it is 
indeed gapless for the pure DM model of spin one. To- 
gether with the presented results, we conclude that the 
quantum phase transition does not occur until D/J goes 
to infinity. 
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FIG. 2: (Color online) The spin gap as the function of D/J. 
The gap function is a smooth function of D/J. Although 
we only compute up to D/J = 2.8, the tendency denies the 
quantum phase transition at finite D/J as the Haldane gap 
case. 

Our result demonstrates a nontrivial consistency be- 
tween the Heisenberg model and the AKLT model. That 
the VBS state in both systems are stable against the DM 
perturbation indicates that they are of the same kind. 
Our results also imply the irrelevance to the Haldane 
gap of the (Si ■ Si+i) 2 term in the AKLT model, since 
the nearest-neighbored DM term can not be gauged away 
due to that term and the spin gap is still robust. 



II. QD MODEL 

In this section, we consider the QJ model given by 
Eq.Q in the square latticed If Q = 0, the ground state 
is the antifcrromagnetic state. If J = 0, the ground state 
is shown to have the VBS order. There is a quantum 
phase transition at J*/Q = 0.04 that both the Neel's 
order and the VBS order vanish, which is shown to be the 
second order phase transition with emergent global U(l) 
symmetry. The small value of the J* /Q indicates that 
the VBS state is very vulnerable when competing with 
the antifcrromagnetic state, that serves an good example 
to explain the current theoretical difficulty that it is very 
difficult to construct a pure local quantum spin model to 
host the VBS state as the ground state. 

Now, we would like to examine if this VBS state is 
also at the weak side when competing with the DM in- 
teraction? In the transition metal oxides, it is quite 
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common that the inversion symmetry is not exactly pre- 
served. Consequently, the DM interaction, regardless big 
or small, is not vanishing. If the VBS state is also very 
vulnerable against the DM interaction, it will be a hard 
task for us to observe it experimentally To answer the 
question, let us take the QJ model plus the DM interac- 
tion and take J = as the example. The Hamiltonian is 
given by the following 



<ijkl> 



Let us call it the QD model and work in the square lattice. 
To perform the calculation, the specific configuration of 
the Dij is needed. Here we borrow the D^ from the 
case of the DM interaction in the perovskite La2CuO j^. 
Although the ground state of La2Cu04 is not the VBS 
state, the configuration of the DM vector in this system 
is quite typical. Readers should be clear that we are not 
studying the system of La2Cu04 here. 
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momentum P. The ground state is always found in the 
5 = and P = sector. However, the first excited state 
is found either in the 5 = and P = sector or in the 
5=1 and P = sector in the calculation for different 
sizes of the system. The first one is referred to the singlet 
excitation, and the latter one is the triplet excitation. We 
expect that the first excitation is due to the broken bond 
of one valence bond and therefore should be in the triplet 
sector. Indeed, we found that the first triplet excitation 
scales faster than the singlet one to have the lower en- 
ergy. Therefore, we look for the first excited state in the 
triplet sector. 
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FIG. 4: (Color online) The finite size scaling of 1/L of the spin 
gap as the function of D/Q. The system size is 4 x L. We 
compute for L = 4, 6, and 8. Both x and the y directions take 
the periodic boundary condition, a is the fitting parameter 
so that the energy gap scales as 1/L a 

Fig. Q is the result of the finite-sized scaling as 1 / L of 
the spin gap for different D/Q values. The system sizes 
are taken to be 4 x L with L = 4, 6, and 8. We fit the 
scaling of the spin gap A(L) by 



FIG. 3: (Color online) The DM vector Dij in the tetragonal 
phase in La2Cu04. Here we only borrow the configuration of 
Dij as a typical example to consider the the DM perturbation 
on the VBS state. 

In the La2CuC>4, the DM interaction occurs in two 
different cases depending on the way of tilting of the CuO 
octahedral. Two situations happen in the orthorhombic 
and the tetragonal phases respectively. Actually, they 
belong to the same kind of configurations, so we only 
study the tetragonal case where the DM vectors D^ are 
given by Fig. 

We apply the Lanczos method to compute the spin 
gap of the Eq.Q. We consider the periodic boundary 
condition in the both x and y directions, and therefore 
the momentum is a good quantum number. The sparse 
Hamiltonian is divided by the subspaces with the quan- 
tum numbers of the total magnetization 5 and the total 



A(L) = A oo +0( — ) 



(9) 



where Aoo is the value by taking L to be infinity. Al- 
though it is not shown here, we also compute the spin gap 
of the pure DM model, namely Q = in the QD model. 
In this case, a = 1, and we found A^Q = 0) = 0. For 
finite Q, a > 1 and change smoothly to a — 1, indi- 
cating the lack of the quantum phase transition at finite 
Q. We also compute the A^ as the function of D/Q 
in the Fig.(j5|. While the error to compute the energy 
using the Lanczos method is given around the order of 
10 _6 (5, the error bars in Fig.fcl are mainly from fitting 
the finite-sized scaling. 

As the DM interaction is turned on, the spin gap re- 
duces. It ramps down rapidly around D/Q ~ 1 but still 
remains finite after reaching the foot of the hill. The re- 
sult can be more accurate if we also do the finite-sized 
scaling along the x direction. Nonetheless, the result will 
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FIG. 5: (Color online) The spin gap Aoo as the function of 
D/Q. The gap does not close at our largest computed D/Q 
= 3.33. 

change quantitatively but not qualitatively. All values 
will become just smaller since the gap due to the finite 
size along the x direction also reduces. Here we compute 
up to D/Q = 3.33 and fail to observe the quantum phase 
transition. From the tendency of the curve, we conclude 
that the gap does not close until Q = 0. 

Let us briefly comment about the introduction of the 
J— term in the QD model, namely the QJD model. Due 
to the reduction of the spin gap by the DM interac- 
tion, we can imagine that the quantum phase transition 
from the VBS state to the AFM state in the QJ model 
should happen at smaller J* in the QJD model. How- 
ever, whether or not the quantum phase transition is still 
of the second order type remains further study. We re- 
mark that it is not an easy task, because the minus-sign 
problem returns in the QJD model, as the QJ model is 
free of that problem. Unfortunately, it is not appropriate 
to apply the Lanczos method to study the QJD model. 



Since it requires a larger system to see the finite-sized 
scaling to locate the transition point, the size that the 
Lanczos method can reach is not enough. 



III. CONCLUSION AND DISCUSSION 

In this paper, we have apply DMRG method and the 
Lanczos method to compute the spin gap in the AKLT 
model and the QD model with the DM interaction. In 
both systems, the spin gap of the VBS state is stable 
against the DM perturbation. The quantum phase tran- 
sition does not occur until the pure DM model is reached. 
Albeit the VBS state is difficult to obtain both experi- 
mentally and theoretically, our results show positively 
that the VBS state will not be destroyed by the DM in- 
teraction which is very common in the transition metal 
oxides or similar chemical families. We understand that 
our result is different from the VBS state of the Heisen- 
berg model in the kagome lattice, where at D/J = 0.1 
the transition from the VBS to the DM ground state oc- 
curs. However, the ground state of the Heisenberg model 
in the kagome lattice has not reach consensus yet. More- 
over, the VBS state in the frustrated lattice could be very 
different from the one in the bipartite lattice. Therefore, 
our results do not cause inconsistency with previous re- 
sults. 

To extend the current research, it is interesting to find 
out the nature of the quantum phase transition from the 
VBS state to the DM ground state. As it is shown in 
the RefpJ, it is the Kosterlitz-Thouless-type transition 
in the Heisenberg model in the spin chain. Although we 
expect the same result in the AKLT model, the one in 
the QD model will be interesting. 
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